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' Aims. We compare six popularly used evolutionary population synthesis (EPS) models through fitting the full optical spectra of six 

representative types of galaxies (star-forming and composite galaxies, Seyfert 2s, LINERs, E+A and early-type galaxies), which are 
taken from the Sloan Digital Sky Survey (SDSS); and we also explore the dependence of stellar population synthesis results on the 
main ingredients of the EPS models; meanwhile we study whether there is an age sequence among these types of galaxies. 
Methods. Throughout our paper, we use the simple stellar populations (SSPs) from each EPS model and the software STARLIGHT 
I to do our fits. Firstly, to explore the dependence of stellar population synthesis on EPS models, we fix the age, metallicity, and initial 

• , mass function (IMF) to construct a standard SSP group. We then use the standard SSP group from each EPS model (BC03, CB07, 

i-^ ' Ma05, GALEV, GRASIL, and Vazdekis/Miles) to fit the spectra of star-forming and E+A galaxies. Secondly, we fix the IMF and 

change the selection of age and metallicity respectively to construct eight more SSP groups. Then we use these SSP groups to fit the 
spectra of star-forming and E+A galaxies to check the effect of age and metallicity on stellar populations. Finally, we also study the 
<^ , effect of stellar evolution track and stellar spectral library on our results. At the same time, the possible age sequence among these 

types of galaxies are suggested. 

Results. Using different EPS models the resulted numerical values of contributed light fractions change obviously, even though the 
dominant populations are consistent. The stellar population synthesis does depend on the selection of age and metallicity, while it 
does not depend on the stellar evolution track much. The importance of young populations decreases from star-forming, composite, 
Seyfert 2, LINER to early-type galaxies, and E+A galaxies lie between composite galaxies and Seyfert 2s in most cases. 
^ ' Conclusions. Different EPS models do derive different stellar populations, so that it is not reasonable to directly compare stellar 

CO I populations estimated from different EPS models. To get reliable results, we should use the same EPS model for the compared 
samples. 
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1 . Introduction mation histories. This is the so-called stellar population synthe- 

sis on galaxies. Two main types of approaches have been de- 
Stellar populations are fundamental characters m revealmg the ei^pirical population synthesis (Faber 1972; Bica 
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formation and evolution of galaxies. The formation of spiral ^^gg. ^^.^^^^ 2OOO; Cid Fernandes et al. 2001) and the 
galaxies IS still hotly debated, and two main channels have been gpg ^ j-^^^^^ j^yg; Bruzual 1983; Worthey 1994; Leitherer 



proposed, either initia collapse of gas at very high redshift in ^ Heckman 1995; Maraston 1998; Vazdekis & Arimoto 1999; 
^ the frame of the tidal torque theory (Eggen, Lynden-Bell & ^^.^^^^j ^ ^.j^^^j^^ 2003; Mai-aston 2005; Cid Fernandes et 
■ ■ Sandage, 1962; White, 1984), or gaseous-nch mergers at inter- 2005). In the empirical population synthesis approach, also 
mediate to high redshifts Hammer et al., 2005, 20(^7, 2009). ^^^^ ' stellar population synthesis with a data base', the ob- 
These two channels for galaxy formation may provide distinct ^^^^^^ ^^^^^^ ^ reproduced by a combination of 
signatures from the analysis of the stellar populations, and it is individual stars or star clusters with different ages and 
relevant to test whether or not stellar population models can be metallicities from a library. The results following this approach 
, '°J^n T- '^^ example Heavens et al., 2004; Panter ^^^^.^^^ ^j^U^^ evolution, and do not allow one to 
et al, 2007). While stars can not be resolved for a majority of ^^j^^ ^^^^^^ ^ ^^^^^j ^ance of galaxies, 
galaxies, therefore many works have been generated on analyz- 
ing stellar populations through the integrated lights of the galax- The EPS approach uses the knowledge of stellar evolution to 
ies. Because the integrated lights hold information about age and model the spectrophotometric properties of stellar populations, 
metallicity disti-ibutions of their stellar populations and stai- for- and has enjoyed more widespread use recently. In this approach, 

the main adjustable parameters are the stellar evolution tracks, 

* email: chenxy@bao.ac.cn the stellar spectral library, the IMF, the star formation history 

** email: ychang@bao.ac.cn (SFH), and the grids of ages and metallicities. EPS is a real phys- 
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ical model, but it is restricted by the lacking of comprehensive 
stellar spectral Ubrary, accurate IMF and SFH, and poor under- 
standing of some advanced phases of stellar evolution, such as 
the blue stragglers (BSs), the horizontal branch (HB) stars, and 
the thermally pulsating asymptotic giant branch (TP- AGE) stars. 

Up to now, several EPS models have been proposed and 
widely used in the stellar population studies on galaxies by ana- 
lyzing their colors, spectra and multi- wavelength spectral energy 
distributions (SEDs), such as BC93 (Bruzual & Chariot 1993), 
Pegase (Fioc & Rocca-Volmerange 1997), GRASIL (Silva et 
al. 1998), GALAXEV (Bruzual & Chariot 2003, BC03), CB07 
(Chariot & Bruzual 2009), SPEED (Jimenez et al. 2004), BaSTI 
(Pietrinferni et al. 2004, Cordier et al. 2007), Ma05 (Maraston 
2005), Starburst 99 (Vazquez & Leitherer 2005), Vazdekis/Miles 
(Sanchez-Blazquez et al. 2006; Vazdekis et al., in preparation), 
GALEV (Anders & Alvensleben 2003; KotuUa et al. 2009), 
SPoT (Raimondo et al. 2005). Some researches have used these 
EPS models to analyze the stellar populations in galaxies and 
star clusters, and even compared them. The results are interest- 
ing, however, most of them focus on analyzing the colors and 
multi-wavelength SEDs of the systems. 

Maraston et al. (2006) used two sets of EPS models to es- 
timate the star formation histories, ages, and masses of seven 
galaxies in the Hubble Ultra Deep Field by analyzing their ob- 
served spitzer mid-IR (the rest-frame-UV) photometry data. One 
of the EPS models is Ma05, which includes the contribution of 
TP-AGB stars, and another one is represented by BC03 (simi- 
lar models are Pegase, Starburst 99 etc.). They concluded when 
they assumed a zero reddening, MaOS gave better fits than BC03 
for these distant passively evolving galaxies at 1.4 < z < 2.7. 
While after dust was included in the fits, Ma05 performed no 
better than BC03. Lee et al. (2007) reconstructed some com- 
posite grids by using BC03, Ma05, and SPoT in the color-color 
diagrams to estimate the average age and metalUcity of spiral 
galaxies. They commented that the scatter among different mod- 
els was large at ages < 2Gyr, and the dominant uncertainties 
arised from the treatment of different evolution phases (e.g. con- 
vective core overshoot, TP-AGB, helium abundance at higher 
metallicities). Longhetti & Saracco (2008) found that the use of 
different models (BC03, Ma05, Pegase, and GRASIL) did not 
bring significant changes in the stellar mass estimates of early- 
type galaxies. Their samples were 10 massive early-type galax- 
ies at 1 < z < 2. Muzzin et al. (2009) fit the UV to NIR SEDs of 
34 K-selected galaxies at z ~ 2.3 (Kriek et al. 2006; 2007; 2008) 
by using BC03, Ma05, and CB07. They concluded that there was 
no significantly better one among them, and the choice of the 
model produced more uncertainties than metallicities. Carter et 
al. (2009) reproduced the optical and near-infrared colours of 14 
nearby elliptical and SO galaxies by using seven different EPS 
models (BC03; Pegase; Starburst 99; GALEV; SPEED; Ma05; 
BaSTI). Their results showed broad agreement on the ages and 
metallicities derived from different EPS models, although there 
were different deviations from the measured broad-band fluxes. 

Conroy et al. (2009a,b,c) carried out a series of works on 
studying the propagation of uncertainties in stellar population 
synthesis modeUng. In their first work they explored the rel- 
evance of uncertainties in stellar evolution, IMF, and stellar 
metallicity distributions to the derived physical properties. They 
subsequently investigated some of the uncertainties associated 
with translating synthetic galaxies into observables, such as stel- 
lar evolution and dust. In their third work they especially per- 
formed the model calibration, comparison, and evaluation, and 
the models included their own flexible stellar population synthe- 
sis model (FSPS), BC03, and Ma05. They found fliat the FSPS 



and BC03 models were able to reproduce the optical and near- 
IR colors of Eh-A galaxies, while the Ma05 model performed 
poorly. They also pointed out significant differences between 
Ma05, BC03 and FSPS in terms of photometry of intermediate- 
age and sub-solar metallicities star clusters (i.e.star clusters in 
the MCs). Namely, their FSPS is better than BC03 and MaOS for 
such cases. 

Similarly, some other works have been applied to star clus- 
ters. Hempel et al. (2005) presented the results of optical and 
near-infrared photometry for globular cluster systems of two gi- 
ant ellipticals. They compared the (V-H) and (V-I) colors from 
BC03 and MaOS, and found the color predictions for a given age 
did not differ significantly, except for metal-poor objects with 
ages < 2Gyr. Pessev et al. (2008) adopted 54 star clusters to eval- 
uate the performance of four EPS models (Vazdekis 1999, BC03, 
MaOS, GALEV) in the optical/near-infrared colour-colour space. 
They argued that each model had strong and weak points and 
there was no model that stood out in all aspects. 

There are also some works on analyzing the spectra of 
galaxies or star clusters. For example, Panter et al. (2007) used 
MOPED to study the star formation history, the stellar mass 
function, and the current stellar mass density in a large sam- 
ple of SDSS DR3 galaxies. They also investigated the effects 
from the choices of the spectral resolution, sky lines, IMF, and 
EPS models (SPEED, Pegase, BC93, MaOS, BC03, CB07) on 
these properties. Their conclusion was that the main impacts on 
estimation of SFH were due to the EPS model, the calibration 
of the observed spectra and the choice of IMF. In 2008, they 
subsequently investigated the impact of model choices by re- 
covering the metallicity history with the same six EPS models, 
and suggested that older EPS models did not produce a clear 
results (Panter et al. 2008). Cid Fernandes and his colleagues 
have worked on spectral synthesis on the SDSS galaxies by using 
the BC03 model and their STARLIGHT code. They also studied 
the properties of the galaxies, such as the dust extinction, stellar 
mass, SFH etc. (Cid Fernandes et al. 2004, 2005, 2007; Mateus 
et al. 2006; Asari et al. 2007, 2009; Stasiriska et al. 2008). 
Koleva et al. (2008) compared different models, and found that 
the Pegase-HR (Le Borgne et al. 2004) and Vazdekis/Miles mod- 
els were in precise agreement, while the BC03 model presented 
biases, which might be due to the poor metallicity coverage of 
STELIB library causing unreliable results at non-solar metal- 
licities. In Koleva et al. (2009a), they used fuU-spectrum fit- 
ting to derive the radial profiles of the SSP-equivalent ages and 
metallicities for a sample of 16 dwarf elliptical galaxies with 
their VLT spectra, and discussed the sensitivity to the population 
model and IMF. Recently, Cid Fernandes & Gonzalez Delgado 
(2009) and Gonzalez Delgado & Cid Fernandes (2009) studied 
the ages, metallicities and dust extinction of 27 star clusters from 
the Magellanic Clouds (from Leonardi & Rose 2003) through 
stellar population synthesis. They adopted STARLIGHT to fit 
their integrated optical spectra in the blue-near-UV range (36S0- 
4600 A). They further compared the combinations of model and 
spectral library. 

Although these impressive progresses in stellar population 
analysis of galaxies and in comparisons among different EPS 
models, we can easily notice that most of them trade on colors 
and multiwavelength SEDs of galaxies. There are no much work 
through fitting the full spectra of galaxies, i.e. the detailed fit- 
tings on their stellar absorptions and continua. Although some 
works do fit spectra to study the properties of galaxies, there are 
no much efforts on comparing the different EPS models through 
fitting spectra on stellar population analysis. Therefore, in this 
work, we will compare six popularly used EPS models carefully 
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through fitting the good quahty optical spectra of a representa- 
tive sample of galaxies taken from SDSS. 

The powerful SDSS has provided good quality spectra of 
hundred thousands galaxies and many other types of astronomi- 
cal objects (such as AGNs) with middle resolution (3 A, /?=2000) 
(Tremonti et al. 2004; Kauffmann et al. 2003; Brinchman et al. 
2004 etc.). These are good samples for reliable studies on stellar 
populations of galaxies from full-spectrum fitting on continua 
and stellar absorption. It also provides a good way to compare 
the different EPS models. 

In this work, we will compare six different EPS models 
(BC03, CB07, Ma05, GALEV, GRASIL and Vazdekis/Miles), 
by doing spectral synthesis analysis on six representative types 
of galaxies from SDSS: 419 star-forming, 326 composite galax- 
ies, 35 Seyfert 2s, 69 LINERs, 502 E-hA galaxies, and 754 early- 
type galaxies. Hence, we check the dependences of the main in- 
gredients of EPS models, i.e. the ages, metallicities, and stellar 
evolution tracks etc., on the stellar population analysis results. 

We will try to answer these questions: Whether these 
different EPS models (with same ingredients) could provide 
same/similar stellar population on one galaxy? If there is dif- 
ference, what is the main ingredient dominates it? What about 
the stellar populations of these representative galaxies? Is there 
any age sequence among them? The latter two questions are 
also very interesting and important, although our main aim is to 
check some aspects of the consistency between models by using 
a single program. 

This paper is organized as follows. In Sect. |2] the detailed 
illustration about the six EPS models and their main ingredients 
are performed. The sample selections are given in Sect. [3] In 
Sect.m we present the methods and the spectral synthesis results 
from different EPS models, and then compare these models. The 
dependences of ages and metallicities on stellar populations of 
galaxies are also checked. In Sect.|5] we discuss the dependences 
of stellar populations of galaxies on stellar evolution tracks, and 
also present the possible age sequence among different types 
of galaxies. Summary and conclusions are given in Sect. |6l 
Throughout this paper we assume the following cosmological 
parameters: Dm - 0.3, Qa - 0.7, and//o=70kms 'Mpc ^ 

2. EPS models 

In this section, we will present the six EPS models which we will 
use and compare in stellar population analysis of galaxies. They 
are BC03, CB07, Ma05, GALEV, GRASIL and Vazdekis/Miles. 
The main ingredients of them contain the stellar evolution tracks, 
the stellar library, the grids of ages and metallicities, the IMF and 
the SFH. 

2.1. The main ingredients of EPS model 

One of the important ingredients in EPS model is the stellar evo- 
lution track, which records the evolution of stars of any mass 
from the zero-age main sequence to later evolutionary stages. 
The tracks with wide range of mass and time are necessary 
for EPS, and many groups have produced public stellar evolu- 
tion models: Padova 1994 (Pa 94, Alongi et al. 1993; Bressan et 
al. 1993; Fagotto et al. 1994 a,b; Girardi et al. 1996) includes 
all phases of stellar evolution from zero-age main sequence 
to the beginning of TP-AGB (for low- and intermediate- mass 
stars) and core-carbon ignition (for massive stars); Padova 1999 
(Pa 99) extends earlier models by the inclusion of the TP-AGB 
phase for stars in the mass range 2Mq < m < 7Mo in ac- 



cordance with the fuel consumption theorem (Schulz 2002); 
Padova 2000 (Pa 00, Girardi et al. 2000); Geneva (Schaller et al. 
1992); Cassisi (Cassisi, Castellani & Castellani 1997; Cassisi, 
DeglTnnocenti & Salaris 1997 and Cassisi et al. 2000); and 
Marigo & Girardi (2007), which is a new synthetic models of 
the TP-AGB evolution. 

The stellar spectral library is another important ingredient 
in EPS model. An ideal stellar library should provide complete 
coverage of the HR diagram, accurate atmospheric parameters 
(Teff, surface gravities log g, metallicities Z etc.), good wave- 
length coverage, good spectral resolution and good calibration. 
Bruzual (2005) and Coelho (2009) reviewed some characteris- 
tics of different stellar spectral libraries. There are two types 
of stellar spectral libraries: empirical library and theoretical li- 
brary. The empirical library (eg. STELIB, Le Borgne et al. 2003; 
UVES POP, Valdes et al. 2004; Indo-US, Jehin et al. 2005; 
MILES, Sanchez-Blazquez et al. 2006; ELODIE, Prugniel et al. 
2007) is based on observations of real stars, so that it is more re- 
liable. However, it is limited by the quality of the observations, 
thus the coverage of parameters are biased towards the typical 
stellar population targeted by the observations. While the theo- 
retical spectral library (eg. Kurucz 1992; BaSeL 1.0, Lejeune et 
al. 1997, 1998) is based on model atmospheres, so that its cov- 
erage of parameters can be at will, but limited by our knowledge 
of the physics of stellar atmospheres. 

The third important ingredient in EPS model is the IMF, 
i.e. the initial distribution of the stars along the main sequence. 
Pioneered by Salpeter (1955), many types of IMF emerged 
(Kennicutt 1983; Scalo 1998; Kroupa 2001; Chabrier 2003; 
etc.). The general form of IMF is <p(M) a M"*'^"'*, and the loga- 
rithmic slope and the upper and lower limits in IMF are particu- 
lar relevant to EPS. 

With the above three main ingredients (stellar evolution 
track, stellar spectral library, and the IMF), we can construct 
SSPs. Then the time-dependent SSPs are convolved with an 
arbitrary SFH to synthesis the spectrum or colors of a galaxy 
(Kennicutt 1998). Throughout this work, we apply a SFH of in- 
stantaneous burst to our analysis i. e. t//(t) = IMq (5(f)- 

2.2. BC03 

A popular library of evolutionary stellar population synthesis 
models, GALAXEV, is computed by using the isochrone syn- 
thesis code BC03 (Bruzual & Chariot 2003). The SSP inod- 
els given by BC03 span a large range of wavelength (91 A ~ 
leOfim, N=6900), age (0 ~ 20Gyr, N=221), and metallic- 
ity Z (0.0001 ~ 0.05, N=6). These models mainly use the 
STELIB/BaSeL3.1 libraries (Le Borgne et al. 2003; Lejeune et 
al. 1997,1998; Westera et al. 2002 and references therein) with 
resolutions of 3A (FWHM) from 3200 A to 9500 A and 20A 
(FWHM) elsewhere. Besides that, BC03 also provides mod- 
els rely on STELIB/Pickles libraries (Pickles 1998). Moreover, 
there are 2 IMFs (Chabrier, Salpeter) and 3 stellar evolution track 
(Pa 94, Pa 00, and Geneva) we can choose. These information are 
given in Table [1] 

2.3. CB07 

CB07 is a new version of BC03, which includes the new stellar 
evolution prescription of Marigo & Girardi (2007) for the TP- 
AGB evolution of low- and intermediate-mass stars. It is given 
by Chariot & Bruzual (2009), and we obtained the model by pri- 
vate communications. The SSP models given by the current ver- 
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Table 1. Comparisons among different EPS models. First, the stellar spectral library, delivering to the spectral resolution, the range 
and number of wavelength, are fixed in our work; Second, in this table, we list the full range of age and metallicity, and we will select 
some of them to construct the groups of SSPs for our fit (Table |2]i; Third, there are several choices of IMFs and stellar evolution 
tracks in some models, and we only list the one we frequently used. 



models 


BC03 


CB07 


Ma05 


GALEV 


GRASIL 


Vazdekis/Miles 


library 


STELIB/BaSeL3.1 


STELIB/Kurucz92 


BaSeL2.0 


BaSeL2.0 


Kuruczl992 


Miles2006 


resolution(A)'' 


3 


3 


20 


20 


20 


2.3 


wavelength (A) 


91-1. 6x10* 


91-3.6x10'* 


91-1.6x10" 


91-1.6x10* 


91-1.2x10' 


3540-7410 


N,, 


6900 


6917 


1221 


1221 


1264 


4300 


age(Gyr)(Number) 


0-20(221) 


0-20(221) 


10-*- 15(67) 


4x10-3-16(4000) 


10-*-20(55) 


0.06-18(50) 


Z(Number) 


0.0001-0.05(6) 


0.0001-0.05(6) 


0.0001-0.07(6) 


0.0004-0.05(5) 


0.0001-0.1(7) 


0.0001-0.03(7) 


IMF 


Salpeter 


Salpeter 


Salpeter 


Salpeter 


Salpeter 


Salpeter 


track 


Pa 94 


Pa94-l-Marigo 07 


Cassisi-l-Geneva 


Pa 99 


Pa* 


Pa 00 



Notes: resolution in visual regions. refer to Bertelli et al. 1994. 



sion of CB07 span a range of wavelength (91 A ~ 36000jum, 
N=6917), age (0 ~ 2QGyr, N=221), and Z(0.0001 ~ 0.05, 
N=6). Moreover, there are two choices of IMF: Chabrier and 
Salpeter These information are given in Table[T] 

2.4. Ma05 

The original version of Ma05 is Ma 98, which is an EPS model 
based on fuel consumption theorem (quite different from other 
models). A grid of SSP models with Z = and an age range 
of 3QMyr ~ ISGyr was constructed. Ma 05 spans a wider range 
of stellar population parameters: 6 metallicities (Z - 0.0001 ~ 
0.07), 67 ages (lOV ~ ISGyr), 2 IMFs (Salpeter, Kroupa), 
and 3 horizontal branch (HB) morphologies (red, intermediate 
or blue; see details in Maraston 2005). However, not each metal- 
licity match all 67 ages: Z = 0.0001 associated to 16 ages 
(1 ~ l5Gyr) with Cassisi tracks; Z = 0.07 associated to 16 ages 
(1 ~ l5Gyr) with Pa 00 tracks; while the rest 4 metallicities as- 
sociated to full 67 ages with Cassisi + Geneva tracks. Moreover, 
the Cassisi tracks are without considering overshooting, while 
the Padova tracks include the effects of overshooting. The stel- 
lar spectra were taken from BaSel library (Lejeune et al. 1998), 
with low spectral resolution, i.e. 5 - lOA up to the visual re- 
gion, 20 - lOOA in the near-IR (wavelength ranges from 91A to 
160;UOT, and A^^i = 1221). These information are given in Table[T] 

2.5. GALEV 

GALEV (GALaxy EVolution) evolutionary synthesis models de- 
scribe the spectral and chemical evolution of galaxies over cos- 
mological timescales, i.e. from the beginning of star formation 
to the present (Kotulla et al. 2009). This code considers both the 
chemical evolution of the gas and the spectral evolution of the 
stellar component, allowing for what they call a chemically con- 
sistent treatment. Thus some SSPs in this model present emis- 
sion lines in their spectra. Additionally the GALEV evolutionary 
synthesis models are interactively available at web-interfacqj 
The SSP models provided by this code cover 5 metallicities 
(0.02 < Z/Zq < 2.5), and 4000 ages (4 x IQ^yr ~ l6Gyr) in 
time-steps of 4Myr. They are based on the spectra from Lejeune 
et al. 1997,1998 (BaSeL 2.0), on 3 IMFs (Salpeter, Scalo, and 
Kroupa), and on the theoretical isochrones from the Pa 99 and 
the Geneva tracks. The lower limit of stellar mass for IMF is 
always Q.IMq, while the upper mass hmits are as follows: for 



Padova, it is set by isochrones (about 5OM0 for super-solar 
metallicity and about IQMq for the rest); for Geneva, it is al- 
ways I2OM0. The range of wavelength is 90A ~ 160// m with 
resolution of 20A in the UV-optical and 50 - lOOA in the NIR 
ranges (Schulz et al. 2002; Anders & Alvensleben 2003). These 
information are given in Table [T] 

2.6. GRASIL 

GRASIL (GRAphite and SILicate) is a population synthesis 
code, which takes into full account the effects of dusty inter- 
stellar medium in galaxy spectra (Silva et al. 1998). It is a 
multi-wavelength model for the combination of stellar popula- 
tion and dust, which absorbs and scatters optical and UV pho- 
tons and emits in the IR-submm region. It is particularly suited 
to the study of the IR properties of dusty galaxies. Additionally, 
GRASIL can be very conveniently run using the web interface 
GALSYNTH 0. The SSP models given by this code spans a 
large range of metallicity (Z = 0.0001 ~ 0.1, N=7), and age 
(lO^yr ~ 20Gyr, N=55). Moreover the models are based on the 
Kurucz (1992) stellar atmosphere model, on 4 IMFs (Salpeter, 
Kennicutt 1983, Miller & Scalo and Scalo) with mass range 
0.15 ~ I2OM0, and on the Padova tracks (Bertelli et al. 1994) 
adding the effects of dusty envelopes around AGB stars (see de- 
tails in Silva's PhD thesis). Besides that, the spectra range from 
91 A to 1200/iOT with resolution of 20 A. These information are 
given in Table [1] 

2.7. Vazdekis/Miles 

The EPS model explored by Vazdekis et al. started from 1996, 
experienced several generations, and predicted for studying old 
and intermediate aged stellar populations. The original version 
V96 (Vazdekis et al. 1996) used the Lick polynomial fitting func- 
tions, which were based on the Lick/IDS stellar library, however, 
the Lick stellar library has not been flux calibrated (V96). Then, 
in 1999 Vazdekis gave an extended version of V96 (Vazdekis 
1999, V99), which provided flux-calibrated spectra. This version 
(V99) predicted SEDs for SSPs in two reduced spectral regions 
in the optical wavelength range (3855 ~ 4476A, 4795 ~ 5465 A) 
at resolution ~ 1.8A (FWHM), with a range of Z (-0.7 < 
logiZ/Ze) < 4-0.2) and ages (1 to 17 Gyr). The input stel- 
lar database is the empirical stellar library of Jones (1999). 
There are 4 types of IMFs provided by this version: Unimodal, 



' http://www.galev.org 



^ |http://galsynth.oapd.inaf.it/galsynth/index.php| 
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Bimodal, Kroupa universal, and Kroupa revised (Vazdekis et al. 
2003 and references therein). Further, a revised version of the 
model of V96 and V99 was explored (Vazdekis et al. 2003), 
which replaced the database of Jones (1999) by Cenarro (2001), 
and replaced the isochrones of Bertelli (1994) by Pa 00. This 
model predicted both the strength of the Ca II triplet feature and 
SEDs in the range 8449 ~ 8952A at resolution 1.5 A (FWHM), 
with metallicities -1.7 < [Fe/H] < +0.2, ages 0.1 ~ l8Gyr, 
and 4 types of IMFs (same as V99). In this code, we adopted 
the latest version based on a new empirical stellar spectral li- 
brary MILES (Sanchez-Blazquez et al. 2006). The model spans 
the range of wavelengths 3540 ~ 7410A at resolution ~ 2.3A 
(FWHM), 7 metallicities -2.32 < log(Z/Zs) < +0.22 with 
Pa 00 tracks, and 50 ages 0.063 ~ 17.78 Gyr (Vazdekis et al., in 
preparation). These information are given in Table[T] 

We summarize the related parameters of all the six EPS mod- 
els in Table[T] which are the most frequent choices in stellar pop- 
ulation analysis and used in this work. 

3. Sample 

We selected six representative types of galaxies as our work- 
ing samples. They are 849 emission-line galaxies including 419 
star-forming galaxies, 326 composite galaxies, 35 Seyfert 2s and 
69 LINERs (taken from Chen et al. 2009); 502 E+A galaxies 
(from Goto 2007); and 754 early-type galaxies (from Hao et 
al, 2006). We downloaded all the ID spectra of them from the 
SDSS, which have been sky-subtracted, then the telluric absorp- 
tion bands were removed before the wavelength and spectropho- 
tometric were calibrated (Stoughton et al. 2002). We corrected 
the foreground Galactic extinction using the reddening maps of 
Schlegel et al. (1998) and then shifted the spectra to the rest 
frame. In order to view the global properties of the galaxies in 
each of the sub-groups, we combined all the spectra in each of 
the sub-groups using the task SCOMBINE in IRAF. That is, we 
combined spectra by interpolating them to a common dispersion 
sampling, and computing the median of the pixels. In this way 
we can also improve well the S/N of the spectra (26 ~ 56 in 
4730A ~ 4780A). Then we work on these six combined spectra. 

3.1. Emission-line galaxies 

The emission-line galaxies were taken from Chen et al. (2009) 
directly. This sample was selected from the main galaxy sam- 
ple of the Sloan Digital Sky Survey (SDSS) DR4 database and 
then was cross-correlated with the IRAS Point Source Catalog 
(PSCz) (5" matching radius). We then selected high signal-to- 
noise ratio (S/N) targets according to the criteria of S/N greater 
than 5cr for H/3,Ha, [W/]/l6583, and S/N greater than 3cr for 
[0III]A5QQ1 . Further the objects with spectra not located at the 
center of galaxies , and those with problematic mask spectra 
were all removed from our sample (see details in Chen et al. 
2009). Finally, 849 objects left in our emission-line galaxies' 
sub-sample. The observed-frame spectral wavelength range is 
3800-9200 A, and the resolution is 3 A( FWHM). 

These emission-line galaxies were further divided into 4 
groups: 419 star-forming galaxies, 326 composite galaxies, 35 
Seyfert 2s and 69 LINERs, according to the emission-line di- 
agnostic diagram (see figure 2 in Chen et al. 2009; Baldwin et 
al. 1981, BPT; Shuder et al. 1981; Veilleux & Osterbrock 1987; 
Kewley et al. 2001; Kauffmann et al. 2003). Composite galaxies 
here refer to objects whose spectra contain significant contri- 
butions from both AGN and star formation (Brinchmann et al. 



2004; Kewley et al. 2006; Chen et al. 2009). Most of them have 
redshift z < 0.1 due to the limit of IRAS observations. 

3.2. E+A galaxies 

E+A galaxies are so called because their spectra look like a su- 
perposition of elliptical galaxies and of A-type stars. On the one 
hand, the features of elliptical are not in point of morphology (al- 
though most of them are elliptical) but in point of non-detection 
of ongoing star formation indications, i.e. these galaxies' spec- 
tra do not show significant emission lines. On the other hand, 
the features of A-type stars are reflected in strong Balmer ab- 
sorption lines. Therefore E+A galaxies have been explained as 
post-starburst galaxies, that is, a galaxy has experienced star- 
burst recently, but truncated it suddenly. Some works have been 
performed on this type of galaxies ( Falkenberg et al. 2009a, b; 
Huang & Gu 2009; Pracy et al. 2009, and references therein). 

We selected our sample of E+A galaxies from the catalogue 
of Goto 2007, which was based on the SDSS DR5, and was an 
extended sample from Goto 2005 on DR2. There were 564 E+A 
galaxies in this catalogue, and the selection algorithm is sum- 
marized as follows: they only used objects classified as galaxies, 
spectroscopically classified not to be a star, and with spectro- 
scopic S/N> 10 per pixel to remove the pollution from nearby 
stars and star-forming regions; then they selected E+A galaxies 
as those satisfied the criteria of H6 equivalent width (EW)> 4AQ, 
[Oil] EW> -2.5A, Ha EW > -3.0A; moreover, they excluded 
galaxies at 0.35 < z < 0.37 due to the sky feature at 5577A. 
From this catalogue, we further discarded 62 E+A galaxies as 
they have problematic mask spectra. Thus our final sample of 
E+A galaxies contains 502 objects. 

3.3. Early type galaxies 

Early-type galaxy has been the most studied objects currently, 
because of its relatively little dust extinction, gaseous interstellar 
medium, and little recent star formation (Vazdekis et al. 1996; 
Cid Fernandes et al. 2009, and references therein). 

Our sample of early-type galaxies was selected from Hao et 
al. (2006), which was based on the SDSS DR4 photometric cat- 
alogue: they selected objects with redshifts smaller than 0.05, 
with the velocity dispersions cover a range of 200 km s"' to 420 
km s ', and not saturated or located at the edge of the corrected 
frame; they further excluded some objects which were unavail- 
able from the SDSS DR4 Data Archive Server (DAS), and those 
with visible dust lanes. Then they visually examined all images 
to confirm the target objects were E/SO galaxies and not con- 
taminated by companion galaxies or bright stars. Moreover, we 
discarded 93 objects with problematic mask spectra, so that our 
final sample of early-type galaxies contains 754 objects. 

4. Comparisons among the six EPS models based 
on spectral synthesis 

4.1. Spectral synthesis method 

We fit the spectral absorptions and continua of the sample 
galaxies to study their stellar populations by using the software 
STARLIGHTQ (Cid Fernandes et al. 2005, 2007; Mateus et al. 
2006; Asari et al. 2007; Chen et al. 2009). It is a program to fit 



^ Absorption lines have a positive EW. 
* http://www.starlight.ufsc.br 



6 X. Y. Chen et al.: Comparing different EPS models through fitting spectra of galaxies 



Table 2. Selections of ages and metallicities for SSPs from each model, which we will use below. BC03 and CB07 are shown in the 
same line. See more details in Sect. 14.21 



models 


ages(Gyr) 


metallicities (Zq) 




Young (< O.lGyr) 


Intermediate (O.lGyr ~ 


2Gyr) 


Old (> 2Gyr) 










Yl 


Y2 


Y3 


Y4 


11 


12 


13 


14 


01 


02 


03 


04 


Zl 


Z2 


Z3 


BC03,CB07 


0.004 


0.010 


0.064 


0.102 


0.286 


0.509 


0.905 


1.434 


3.00 


6.00 


10.00 


13.00 


0.2 


1.0 


2.5 


Ma05 


0.004 


0.010 


0.065 


0.100 


0.300 


0.500 


0.900 


1.500 


3.00 


6.00 


10.00 


13.00 


0.5 


1.0 


2.0 


GALEV 


0.004 


0.012 


0.064 


0.100 


0.280 


0.500 


0.900 


1.432 


3.00 


6.00 


10.00 


13.00 


0.2 


1.0 


2.5 


GRASIL 


0.004 


0.010 


0.060 


0.100 


0.300 


0.500 


0.900 


1.500 


3.00 


6.00 


10.00 


13.00 


0.2 


1.0 


2.5 


Vazdekis/Miles 






0.063 


0.100 


0.280 


0.500 


0.890 


1.410 


3.16 


6.31 


10.00 


12.59 


0.2 


1.0 


1.5 



an observed spectrum 0,i with a model that adds up to Nt 
SSPs with different ages and metallicities from different stellar 
population synthesis models. A Gaussian distribution centered 
at velocity v, and broadened by cr^ models the line-of-sight stel- 
lar motions. The fit is carried out with the Metropolis scheme 
(Cid Fernandes et al. 2001), which searches for the minimum 
= - M^)aj^]^, where the reciprocal of weight w^' is 

the error in Oa except for masked regions. Pixels that are more 
than 3cr away from the rms Oa-M^ are given zero weight by the 
parameter 'clip'. The STARLIGHT group has carefully checked 
the reliability of this software by analyzing the "stellar popula- 
tions" of fake galaxies made by known SSPs (see figure4 in Cid 
Fernandes et al. 2005, and figurel in Cid Fernandes et al. 2004). 

Throughout our fit, we used the reddening law of Calzetti et 
al. (1994, CAL hereafter). The fitted wavelength range was from 
3700 to 7400 A. A power law stands for non-stellar component 
(Koski et al. 1978) Fy a y^'^ was added when we fit the spectra 
of Seyfert 2s and LINERs. The mask regions were the same as 
Chen et al. (2009), which include the strong emission-lines and 
the bad pixels. Additionally, we do not consider different weights 
in different regions of the spectra, since we have tested this and 
found that it did not affect our spectral synthesis results. 

In the outputs of STARLIGHT, one of the most important pa- 
rameters to present stellar population is the population vector x. 
The component Xj(y = 1, ...,Nt) represents the fractional contri- 
bution of the SSP with age tj and metallicity Zj to the model flux 
at the normalization wavelength Aq - 4020A. Equivalently, an- 
other important parameter, the mass fraction fij, has the similar 
meaning. 

4.2. SSP selections 

As mentioned in Sect.|2] we select 12 (or 10) representative ages 
and 3 metallicities from each EPS model to construct different 
SSP groups to do our fits, and we list the details in Table |2] 
We arrange the ages of SSPs into 3 bins: young populations 
with age < O.lGyr, intermediate-age populations with age be- 
tween 0.2Gyr and 2Gyr, and old ones with age > 2Gyr. We 
note that these criteria are different from Chen et al. (2009), 
which followed the work of Kong et al. (2003) (i.e. young with 
age < O.SSGyr, old with age > lOGyr, and intermediate with 
ages between these two). Since the width of Balmer absorption 
line can indicate the young population with age < IGyr, while 
the ones older than 1 .5 - 2Gyr will be dominated by emission 
lines. Generally from SDSS spectra it is uneasy to disentangle 
3 - 5Gyr population from 10G3'r populations, and Mathis et al. 
(2006) also pointed out that the signatures of intermediate-age 
stars (0.5 ~ 4Gyr) are masked by those of younger and older 
stars, so we adopt the new criteria here. 

The logic of our work is in five steps as follows: 



1. First, we fix some main ingredients for the EPS models, i.e., 
the stellar library, the stellar evolution track and the IMF 
(Salpeter) are selected according to Table [1] An instanta- 
neous burst SFH is always applied. 

2. Then we select the first SSP group including 6 SSPs at 6 dif- 
ferent ages and at solar metallicity (see Table |2]i from each 
EPS model. It is called group No.(l) and as our "standard 
case" (see Table [3]i: two young SSPs (¥1=0.004, Y2~0.01 
Gyr, except for Vazdekis/Miles, in which the youngest ages 
Y3~0.06 and Y4~0.10 Gyr are selected), two intermedi- 
ate SSPs (12-0.5, 13-0.9 Gyr) and two old SSPs (O3~10, 
04~13 Gyr). This "standard case" will be used to analyze 
the stellar populations of star-forming and E-i-A galaxies in 
Sect. 1431 

3. Next, we construct another three SSP groups by 
adding/removing some young, intermediate and old 
SSPs at solar metallicity. Then we can check the age 
dependences. The corresponding three groups are: 
No.(2): 3 SSPs with ages of Yl (Y3 for Vazdekis/Miles), 12, 
and 03 at Zq (Z2); 

No.(3): 9 (or 8) SSPs with ages of Yl, Y2, Y3, 12, 13, 14, 
02, 03, 04 (only Y3 and Y4 for the young populations in 
Vazdekis/Miles); 

No.(4): 12 (or 10) SSPs with ages of Yl, Y2, Y3, Y4, II, 
12, 13, 14, Ol, 02, 03, 04 (only Y3 and Y4 for the young 
populations in Vazdekis/Miles); 

4. After that, we construct another three SSP groups by adding 
SSPs with sub-solar metallicity or/and super-solar metallic- 
ity based on group No.(l). Then we can check the metallicity 
dependences. These three groups are: 

No.(5): 12 SSPs with 6 ages in common with the group 
No.(l) but at two metallicities: Zl (0.2 or 0.5 Zq) and Z2 

(Zo); 

No.(6): 12 SSPs with 6 ages as in group No.(l) but at two 
metallicities: Z2 (Z©) and Z3 (2.0, 2.5 or 1.5 Zq); 
No.(7): 18 SSPs with 6 ages as in group No.(l) but at three 
metallicities: Zl, Z2 and Z3. 

5. Finally, we fix the ages and only change the metallicity of 
group No.(l) to build the last two SSP groups. Such that we 
can check the age-metallicity degeneracy in simple way. 
No.(8): 6 SSPs with 6 ages as in group No.(l) but at sub- 
solar metallicity (Zl); 

No.(9): 6 SSPs with 6 ages as in group No.(l) but at super- 
solar metallicity (Z3). 

Table|3]list the details about the 9 SSP groups No.(l)-(9). 

4.3. Spectral analysis and comparisons among six different 
EPS models 

Among the six EPS models, the SSPs in three of them (BC03, 
CB07, Vazdekis/Miles) have similar spectral resolution, and are 
comparable to the resolution of SDSS spectra; another three 
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Table 3. Seven SSP groups used in our stellar population analysis. Please see Table|2]for the meanings of the symbols. "VM" refers 
to the model Vazdekis/Miles. Group No.(l) is our "standard case", and is the basis of other groups. 



Models 


SSPs 


Ages 


Metallicities 


Notes 


No.(l) 


6 SSPs 


Yl Y2; 12 13; 03 04 
Y3 Y4; 12 13; 03 04 


Z2 
Z2 


for VM 


No. (2) 


3 SSPs 


Y1;I2; 03 
Y3; 12; 03 


Z2 
Z2 


for VM 


No.(3) 


9 SSPs 
8 SSPs 


Yl Y2 Y3; 12 13 14; 02 03 04 
Y3 Y4; 12 13 14; 02 03 04 


Z2 
Z2 


for VM 


No. (4) 


12 SSPs 
10 SSPs 


Yl Y2 Y3 Y4; 11 12 13 14; 01 02 03 04 
Y3 Y4; 11 12 13 14; Ol 02 03 04 


Z2 
Z2 


for VM 


No.(5) 


12 (6x 2) SSPs 


Yl Y2; 12 13; 03 04 
Y3 Y4; 12 13; 03 04 


Z2Z1 
Z2 Zl 


for VM 


No. (6) 


12 (6x 2) SSPs 


Yl Y2; 12 13; 03 04 
Y3 Y4; 12 13; 03 04 


Z2 Z3 
Z2 Z3 


for VM 


No. (7) 


18 (6x 3) SSPs 


Yl Y2; 12 13; 03 04 
Y3 Y4; 12 13; 03 04 


Zl Z2 Z3 
Zl Z2 Z3 


for VM 


No.(8) 


6 SSPs 


Yl Y2; 12 13; 03 04 
Y3 Y4; 12 13; 03 04 


Zl 
Zl 


for VM 


No. (9) 


6 SSPs 


Yl Y2; 12 13; 03 04 
Y3 Y4; 12 13; 03 04 


Z3 
Z3 


for VM 



(Ma05, GALEV, GRASIL) provide lower spectral resolution 
SSPs, ~20A. Thus, when we use SSPs in the latter three EPS 
models, we need to decrease the resolutions of the observed 
spectra to match the models. We have tried three methods to do 
that, the first is the GAUSS task in IRAF, the second is DisgallD 
developed by our French collaborators, and the third is re-bin 
program kindly provided by Dr. Jingkun Zhao at NAOC. We 
find there are no much difference among the results from these 
three methods. Therefore, finally we adopt the GAUSS task in 
IRAF for de-resolution on our spectra. 

In this section, we adopt SSP group No.(l) from 6 different 
EPS models to analysis the stellar populations. In that way, the 
selection of age, metallicity, and IMF are fixed, and the only free 
factor is EPS models. So that we can check the dependences of 
stellar population synthesis on EPS models. Take star-forming 
and E-hA galaxies as examples, we analyze their stellar popula- 
tions by using group No.(l) from these six EPS models, and then 
compare the results. 

4.3.1. Star-forming galaxies 

In Fig [1] we present the spectral fitting results and the light frac- 
tion distributions for the combined spectra of star-forming galax- 
ies by using six SSPs (group No.(l) in Table|3]l from six differ- 
ent EPS models. Every three panels as a group correspond to 
the results from one EPS model. The top panel shows the ob- 
served spectra (black line), the synthesis spectra (red line), and 
the error spectra (green line); the middle one shows the resid- 
ual spectra; and the bottom one shows the light fraction dis- 
tribution. From the top two panels in each model, we can see 
that the synthesis spectra match well with the observed spec- 
tra in most cases. However, one point we would like to notice 
is, there is a trough around HjS when we use BC03 and CB07, 
and it disappears when we use Ma05, GALEV, GRASIL and 
Vazdekis/Miles. Asari et al. (2007) suggested that this could re- 
late to calibrations in the STELIB library in this spectral range, 
and we also confirmed it in our earlier work (Chen et al. 2009) 
by using star clusters. We speculate that the bad fit may also be 
due to the way to handle the shape of the continuum in Starlight. 
To check this point, we further tried another program, ULySS 
(Koleva et al. 2009b), to perform similar fits, and we got much 



better fits around Hfi. However, this problem does not affect 
much on our stellar population analysis. Besides that, we find 
the Ma05, GALEV and GRASIL models seem to be superior 
in reproducing the important Balmer lines and continuum. This 
could be due to the decreasing resolution, which reduces the de- 
tails in spectra. 

The bottom panels in Fig[T]show that the resulted stellar pop- 
ulations from different EPS models are not exactly same for the 
star-forming galaxies, although the dominant populations are all 
young plus intermediate populations. 

The top line of Table |4] presents the contributed light frac- 
tions from the young (Y), intermediate (I) and old (O) popu- 
lations. The related numbers reveal that star-forming galaxies 
are composite of young, intermediate and old populations. The 
differences of the contributed light fractions among these EPS 
models are in a range of 0.01% ~ 18.75%. Vazdekis/Miles is 
different from others, and it results in 58% (27% more than that 
of CB07) young and 20% intermediate populations. The reason 
of this could be Vazdekis/Miles has no enough young popula- 
tions as others, and its youngest age is ¥3=0.063 Gyr, older than 
others, which may result in higher fraction of young population. 
We believe that the error-bars of these resulted light fractions 
are small. As mentioned by Cid Fernandes et al. (2005) that they 
presented the error-bars centered on the mean values obtained by 
fitting 20 realizations of each of 65 test galaxies. And their three 
condensed populations are well recovered by Starlight, with un- 
certainties smaller than 0.05 (young: t < 10**), 0.1 (intermediate: 
10'^ < t < \Q\ and 0.1 (old: t > 10'') for S/N > 10 respectively. 
In our fittings, the code provides the values of last-chain-values 
for 7 chains, and we find that most of the discrepancies in these 
adopted values are less than 1%. So we will not consider error- 
bars in our following studies. 

Moreover, it also shows that BC03 and CB07 result in -30% 
young stellar population and ~50% intermediate age one. Ma05, 
GALEV and GRASIL result in comparable young and inter- 
mediate populations (both ~44% in Ma05) or slightly higher 
young populations than the intermediate ones (-47% vs. 39% in 
GALEV; -50% vs. 36% in GRASIL). In other words, there is an 
inversion of the dominant population between BC03, CB07 on 
one side, and GALEV, GRASIL and Vazdekis/Miles on the other 
side. In order to explain this phenomenon, we directly compare 



8 



X. Y. Chen et al.: Comparing different EPS models through fitting spectra of galaxies 



wavelength (A) 
5000 6000 



wavelength (A) 
5000 6000 



wavelength (A) 
5000 6000 



'(a) BC03 


' ■ 1 ' 






,1 1 1 1 1 1 1 1 1 1 1 1 1 1 

: 1 1 1 1 , 1 1 , , , 1 , 1 , ; 


' 1 ' 
1 


' 1 ' 
1 





(b) CB07 






' ' 1 ' . 






: , 1 , , , , 1 , , , , 1 , , , , 1 , -i 






J , ^ 



log age fyr] 



wavelength (A) 
5000 6000 



log age [yr] 

wavelength (A) 
5000 ■ 6000 



log age [yr] 

wavelength (A) 
lOOO ■ 6000 




og age [yr] 



Fig. 1. Spectral fitting results of star-forming galaxies by using six SSPs (group No.(l)) from different models: (a) BC03, (b) CB07, 
(c) Ma05, (d) GALEV, (e) GRASIL and (f) Vazdekis/Miles. In each EPS model, top panels: comparison of synthesis spectrum (red 
line) with the observed spectrum (black line), and green line shows the error spectrum; middle panels: the residual spectrum; bottom 
panels: distributions of light fractions. 



Table 4. Stellar populations in star-forming and E+A galaxies by using different EPS models with SSP group No.(l) (Tables [2] EJl. 
Y represents young populations, I represents intermediate populations, O represents old populations. Please see details of age and 
Z selections for each SSP group in text. 



Types 


SSP groups 


Age bin 


BC03 


CB07 


Ma05 


GALEV 


GRASIL 


Vazdekis/Miles 


star- 


No.(l) 


Y 


35 


31 


44 


47 


50 


58 


forming 


6 SSPs 


I 


54 


51 


44 


39 


36 


20 


galaxies 


Z=Z2 


O 


11 


18 


12 


14 


14 


22 


E+A 


No.(l) 


Y 


13 


9 


13 


22 


20 


19 


galaxies 


6 SSPs 


I 


77 


75 


79 


71 


70 


65 




Z=Z2 


O 


10 


15 


8 


7 


10 


16 



the 6 SSPs from these models, which we used in the fittings. 
And we find the young SSPs of BC03 and CB07 hold relatively 
higher fluxes at blue band than others when we normalize them 
at around 6000A. Take Y2 as an example, GALEV and GRASIL 
as one group, the fluxes at blue band of which are lower than that 
of BC03 and CB07, and Ma05 is between these two groups. This 
explains the inversion of young and intermediate populations be- 
tween the two groups, and further Ma05 produces comparable 
young and intermediate ones. As for Vazdekis/Miles, the two 
SSPs in young age-bin are older than others, thus the resulted 
young population is much higher than the intermediate one. The 
reason to cause such difference in SSPs among these EPS models 
could be their different stellar libraries. STELIB is used by BC03 
and CB07, but BaSeL2.0 or Kurucz are used in GALEV, Ma05 



and GRASIL. The former is empirical from real observed stars' 
spectra with higher resolutions and the latter are theoretical ones 
with lower resolutions. The limit of STELIB library has been 
fully discussed by Koleva et al. (2008) and Gonzalez Delgado 
& Cid Fernandes (2009). Another possible reason causing such 
difference between the two groups could be the spectral resolu- 
tion. Since BC03 and CB07 have higher resolutions (3A), while 
Ma05, GALEV and GRASIL have lower resolutions (20A). To 
check this factor, we degraded the 6 SSPs in BC03 from 3Ato 
20A, and then performed similar fits. We find that the discrepan- 
cies between the two groups decrease then. 
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Fig. 2. Same as Fig. [T] but for E+A galaxies. 
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4.3.2. E+A galaxies 

In Fig. |2] we perform the spectral fitting results and the light 
fraction distributions of E+A galaxies by using six SSPs (group 
No.(l) in Table O from six different EPS models. 

We find that the matching between the synthesis spectra and 
the observed spectra is good. Meanwhile, the troughs around 
//g become shallower associated to BC03 and CB07 than that 
of star-forming galaxies, and disappear in panels corresponding 
to Ma05, GALEV, GRASIL and Vazdekis/Miles, which further 
confirm the conclusion we draw in studying star-forming galax- 
ies. Similarly, the bottom panels show that different EPS models 
do not give exactly same populations for E+A galaxies, even 
though the peaks of distributions are more concentrated. 

In the bottom line of Table |4] we also list the resulted 
light fractions of E+A galaxies. It shows that the contributions 
from different models vary within a range of 0.14% ~ 13.53%. 
Moreover, the percentages of intermediate-age populations are 
all greater than 70% except for Vazdekis/Miles (65.45%). Thus 
we suggest that E+A galaxies are dominated by intermediate- 
age populations, and are composed of young, intermediate-age 
and old populations. 

In all, to sum up the points about stellar population analy- 
sis among different EPS models which we have just indicated, 
among different EPS models the numerical values of light frac- 
tions change obviously, although the dominant populations are 
consistent. Additionally, it seems that EPS models still have dif- 
ficulties in disentangling stellar populations, i. e. for a given 
galaxy, different models would give different order of impor- 
tance of old, intermediate and young stellar population (espe- 



cially in the case of star-forming galaxies). We suggest that 
the use of near-IR photometry maybe a way to improve this. 
However, as mentioned by Eminian et al. (2008), Maraston 
(2005) and Conroy et al. (2009a), TP-AGB stars dominate the 
near-IR light of galaxies. So we should be careful that the uncer- 
tainties associated with the TP-AGB phase will probably hamper 
significant improvements. 

4.4. Dependences on ages and metallicities 

In this section we will investigate the dependence of stellar pop- 
ulation synthesis on the selection of age and metallicity. As men- 
tioned in Sect. 14.2! Table [2] and Table [3] we adopt nine SSP 
groups in total, No.(l)-(9). 

4.4.1. Age dependences - SSP groups No.(2)-(4) in Table|3] 

As mentioned above, we fix the metallicity as solar one, then 
we only change the selections of age to construct SSP groups 
No.(2)-(4). Then we use SSP groups No.(2)-(4) from different 
EPS models to analyze the stellar populations of star-forming 
and E+A galaxies. The resulted stellar populations are compared 
with each other and also with group No.(l), so that we can check 
the age effects on stellar population analysis. 

The top panels in Fig. [3] show the light fractions of young 
(top-left panel), intermediate (top-middle panel), and old (top- 
right panel) populations for star-forming galaxies by using the 6 
EPS models. In each panel, the horizontal axis represents the 6 
EPS models, and the vertical axis represents the contributed light 
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Fig. 3. The light fractions' contributions to star-forming galaxies in bins of young population (left three panels), intermediate-age 
population (middle three panels), and old population (right three panels). The x-axis represents six EPS models, and the y-axis 
represents the light fractions in percentage. Different lines connect points from different SSP groups. Top panels: solid line stands 
for group No.(l), dotted line stands for group No. (2), short-dashed line stands for group No. (3), and long-dashed line stands for 
group No. (4). Middle-panels: solid line stands for group No.(l), dotted line stands for group No. (5), short-dashed line stands for 
group No. (6), and long-dashed line stands for group No. (7). Bottom-panels: solid line stands for group No.(l), dotted line stands 
for group No. (8), short-dashed line stands for group No. (9). 
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Fig. 4. Same as Fig.|3] but for E-i-A galaxies. 



fractions. The solid lines connect the resulted light fractions cor- 
responding to group No.(l). The dotted, dashed and long dashed 
lines refer to the connections of the resulted light fractions from 
the SSP groups No. (2), (3) and (4), respectively. 

From these comparisons, generally we find that adding more 
SSPs within each age bin (based on the "standard case") results 
in more young populations(see the top-left panel), but remov- 
ing SSPs does not result in consistent trend for them. In each 
EPS model, most of the changes in young populations (top-left 
panel) are ~1-10% among different SSP groups. The biggest 
changes occur between the groups No.(l) and No. (4) in CB07 
(~16%) as well as between groups No. (2) and No. (4) in GALEV 
(~15%). As to the general trend of intermediate-age popula- 
tions (the top-middle panel), adding/removing SSPs results in 
less contributions. The most obvious changes occur between the 
groups No.(l) and No.(4) in CB07 (-20%) and between No.(2) 
and No. (3) in Ma05 (~ 14%). The changes in old populations (the 
top-right panel) are related to the young and intermediate ones, 
since the sum of these three populations is equal to 100%. 

Similarly, we plot the corresponding results for E-i-A galax- 
ies in top panels of Fig. |4]to check the age dependences. The 



symbols and lines have the same meanings as in Fig. [3] It con- 
firms the general trend in star-forming galaxies, i.e. adding more 
SSPs increases the young populations and decreases the inter- 
mediate populations. And 3 SSPs result in complex changes for 
the different models. In details, the top-left panel shows that 
within each EPS model, most of the changes in young popula- 
tions are around 0-5% among different SSP groups. The biggest 
ones occur between groups No.(l) and No. (4) in Ma05 as well 
as between groups No. (2) and No. (3) in GALEV (~ 20 %); the 
changes in intermediate-age and old populations are obvious in 
CB07 and Ma05, up to ~ 18% and ~ 25%. The reason could 
be that CB07 and Ma05 specially include the TP-AGB contribu- 
tions to build their SSPs, and the E-i-A galaxies are right domi- 
nated by intermediate-age populations, which are quite related to 
the TP-AGB stars. This could also explain the most obvious dis- 
crepancies in the intermediate-age populations of star-forming 
galaxies with different models (see top-middle panel of Fig.[3]l. 
Moreover, the obvious discrepancy of group No. (2) with 3 SSPs 
from others may mean that 3 SSPs are not enough for such stellar 
population analysis on galaxies. 
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In summary, from the top panels in Figs. [3] HI we can com- 
ment that the SSPs with different sets of ages present obviously 
different results (such as CB07, Ma05), although the discrepan- 
cies from both the models and ages will not change the dominant 
population of the galaxies. 

4.4.2. Metallicity dependences - SSP groups No.(5)-(9) in 
Table |3] 

We fix the selections of age as the same as in the group No.(l) 
(Table |3] solar metallicity), then we add sub-solar metallic- 
ity or/and super-solar metallicity to construct groups No.(5)-(7) 
(Table|3]l. By comparing the results from these SSP groups with 
those from group No.(l), we can check the metallicity effects on 
spectral synthesis results. 

The middle panels in Fig. |3] show the light fractions' con- 
tributions from young (middle-left panel), intermediate (middle- 
middle panel), and old (middle-right panel) populations to star- 
forming galaxies by using the 6 EPS models. The solid lines con- 
nect the resulted light fractions corresponding to group No.(l) 
from different EPS models. The dotted, short dashed and long 
dashed lines refer to the connections of the resulted light frac- 
tions from the SSP groups No. (5), (6) and (7), respectively. 

These results show that within each EPS model all of the 
changes in young populations are less than 5%. The changes in 
intermediate-age populations (middle-middle panel) are larger 
and are in a range of 0-17%, such as ~ 17% between groups 
No.(5) and No.(6) in CB07, ~ 13% between groups No.(5) and 
No. (6) in Ma05. For the old population (middle-right panel), the 
changes are also obvious among these SSP groups (0-17%), e.g., 
~ 14% between groups No.(5) and No.(6) in CB07, ~ 17% be- 
tween groups No. (5) and No. (6) in GALEV. Similarly, the more 
obvious changes in CB07 and Ma05 may be related to the in- 
volved TP-AGB stars in their SSPs. 

Similarly, we present the corresponding results for E+A 
galaxies in the middle panels in Fig. H) It shows that within 
each EPS model, all changes in young populations are not ob- 
vious (less than 6%); generally the changes in intermediate-age 
populations and old populations are also not obvious except for 
CB07. The differences is about 15% between groups No. (5) and 
(6) in intermediate population, and about 15% between groups 
No.(5) and (7) (as well (6)) in old population (within CB07). 
These could also be due to the TP-AGB effect. 

In a word, from the middle panels in Figs. [3] ID we suggest 
that SSPs with different selections of metallicities give different 
results. 

By comparing the top panels with the middle panels in 
Figs. [3] m We can see that the effects of metallicity is less than 
ages. However, this phenomenon may be due to either the way 
we show our results or the age-metallicity degeneracy. Note that 
up to now, we focused on the contributions of light fractions. 
While other relevant aspects, such as M/L, maybe sensitive to 
metallicity. Maraston (2005) has shown some parameters' de- 
pendences on metallicity, for example: the IR indices C2 in their 
figure 16, the M* /L in their figure 23, the D4000 in their fig- 
ure 25, and the indices CaT*, CaT and PaT in their figure 26. 
Longhetti & Saracco (2009) also performed the M/L as a func- 
tion of metallicity in their figures 3 ~ 5. 

Another notice is that there could be age-metallicity degen- 
eracy following our adding or removing SSPs in the top and mid- 
dle panels in Figs. |3] |4] For this reason we constructed the last 
two SSP groups, in which we only changed the metallicities and 
fixed all the ages. We present the results of these two groups 
in the bottom three panels in Figs. [3] H) From the bottom three 



panels we see a much larger scatter (up to 40%) among differ- 
ent SSP groups than that of middle panels. Such that we have to 
point out the age-metallicity degeneracy exists in our Figs. [3] E] 
In our current paper, we will not do more work on this classic 
degeneracy between age and metallicity, and some related stud- 
ies can be found in Cid Fernandes et al. (2005), Eminian et al. 

(2008) , Carter et al. (2009), Cid Fernandes & Gonzalez Delgado 

(2009) . 

5. Dependences on evolution track, stellar spectral 
library, and the age sequence of galaxies 

In this section, we further check the effects of stellar evolution 
tracks on the stellar population analysis of galaxies. We take 
BC03 and GALEV models as examples to do such check. We 
will use these two models to analyze the stellar populations of all 
the six representative samples of galaxies, i.e. the star-forming, 
composite galaxies, Seyfert 2s, LINERs, E-hA and early-type 
galaxies. Hence we can compare the different stellar populations 
in each galaxy type. 

5. 1 . Dependences on stellar evolution track and stellar 
spectral library 

We use SSP group No.(l) (in Table [3]l with Salpeter IMF, and 
with two different sets of stellar evolution tracks (Pa94 and PaOO 
in BC03; Pa99 and Geneva in GALEV). The contributed light 
fractions of young, intermediate and old populations are listed 
in Table |5] The first two columns corresponding to BC03 shows 
that within each galaxy type, among different choices of stellar 
evolution tracks, the corresponding changes in young popula- 
tions are 0-4%, in intermediate populations are 0-9%, and in old 
populations are 0-10%. While most of the results of GALEV 
(last two columns) are similar to BC03, except for the changes 
in intermediate and in old populations of early type galaxies 
( 20%). These results imply that changing stellar evolution track 
will not change much the resulted stellar populations of all the 
six types of galaxies. 

Additionally, we mentioned in Sect.4.3.1 that different se- 
lections of libraries may be another factor effects our results. 
Therefore, we compared the two different libraries provided by 
BC03 (Table|5j. One library is STELIB and another is BaSeL3. 1 , 
and we fix the IMF as Salpeter and the track as Pa94. We find that 
the results of BaSeL3. 1 is different from that of STELIB, but the 
results of BaSeL3. 1 is similar with that of GRASIL and GALEV. 
This supports our previous discussions in Sect.4.3. 1 . Besides our 
conclusions, Gonzalez Delgado & Cid Fernandes (2009) com- 
mented that the STELIB-based models produced metallicities of 
the star clusters systematically smaller by about 0.6 dex with re- 
spect to that was found with other models, but ages were proba- 
bly right. The fact was that STELIB significantly underestimated 
metallicities with respect to both the Call triplet and the metal- 
licities from Leonardi & Rose (2003). A similar conclusion was 
derived by Koleva et al. (2008) in their analysis of Galactic glob- 
ular clusters. 

5.2. The age sequence 

Since all the six types of sample galaxies are studied on their 
stellar populations, we can further study whether there is any 
age sequence among them. This could be an extended study of 
Chen et al. (2009), in which we have shown that there is an age 
sequence from star-forming galaxies, through composite galax- 
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Table 5. Stellar populations in each galaxy type by using six 
SSPs (group No.(l)) from BC03 (columns 3-4) and GALEV 
(columns 5-6) in different cases of tracks, and the results from 
BaSeL3.1 of BC03 (column 7). PL represents power law, and 
other symbols are the same as Table |4] 



galaxy types 


age 


BC03 


GALEV 






bin 


Pa94 


PaOO 


Pa94 


Geneva 


BaSeL3.1 


star-forming 


Y 


35 


31 


47 


41 


49 




I 


54 


51 


39 


48 


40 




O 


11 


18 


14 


11 


11 


composite 


Y 


31 


29 


42 


39 


40 




I 


53 


50 


39 


47 


42 




O 


16 


21 


19 


14 


17 


E+A 


Y 


13 


9 


22 


20 


19 




I 


77 


71 


71 


72 


72 




O 


10 


20 


7 


8 


9 


Seyfert 2 


Y 


6 


10 


2 










I 


32 


23 


34 


30 


29 




O 


35 


42 


26 


31 


29 




PL 


27 


25 


38 


39 


42 


LINER 


Y 





3 













I 


24 


16 


26 


22 


20 




O 


58 


65 


50 


58 


54 




PL 


18 


16 


24 


20 


26 


early type 


Y 


6 


7 













I 


7 


2 


40 


20 


31 




O 


87 


91 


60 


80 


69 



ies, Seyfert 2s to LINERs, since the young populations decrease 
following this sequence. 

In addition to the information about the dependences of spec- 
tral synthesis results on stellar evolution track, the numbers in 
Table |5] also shows us a possible age sequence among six dif- 
ferent galaxy types. In our earlier work, we have suggested that 
there was an age sequence: from star-forming, composite galax- 
ies, Seyfert 2s to LINERs. Here we confirm this results, and we 
add E-i-A and early-type galaxies. It shows that the young popu- 
lations decrease from star-forming, composite galaxies, through 
E-hA galaxies, Seyfert 2s, LINERs to early-type galaxies. And 
E-hA galaxies are dominated by intermediate-age populations, 
and the intermediate-age populations in E-i-A galaxies are far 
more than others. While early-type galaxies are certainly the old- 
est one among all galaxy types. Therefore, we conclude there 
is a possible age sequence from star-forming, composite galax- 
ies, through E+A galaxies, Seyfert 2s, LINERs to early-type 
galaxies. A similar conclusion was reached by Schawinski et al. 
2007, who identified an evolutionary sequence from star forma- 
tion via nuclear activity to quiescence, i. e. from star-forming 
via transition region (something like the composite galaxies in 
our sample) and Seyfert AGN and LINER to quiescence. This is 
also consistent with Cid Fernandes et al. (2009b) and Stasinska 
(2008). 

6. Conclusion 

We compare 6 popularly used EPS models (BC03, CB07, Ma05, 
GALEV, GRASIL, Vazdekis/Miles) in this work. We use the 
SSPs they provide to fit the full optical spectra of six represen- 
tative types of galaxies (star-forming and composite galaxies, 
Seyfert 2s, LINERs, E+A and early-type galaxies), which are 
taken from SDSS. In details, we investigate the dependences of 
stellar population synthesis on EPS models, age, metallicity, and 
stellar evolution track. We also study the age sequence of these 
different types of galaxies. 



1 . The spectral fits of various galaxies by using different EPS 
models are excellent except possibly problematic regions, 
such as the area around Hp line in few cases. 

2. We select 6 SSPs with fixed ages, metallicities, IMF and stel- 
lar evolution tracks from each EPS model to fit the spectra of 
star-forming and E+A galaxies. So that we can explore the 
dependences of stellar population synthesis on EPS models. 
We comment that it is a complex job to do stellar popula- 
tion analysis on galaxies, different EPS models may result in 
quite different results (such as BC03 vs. Vazdekis/Miles). 

3. We fix the IMF, metallicity, stellar evolution tracks, and 
change the selection of ages to construct different SSP 
groups from the 6 EPS models, and then fit the spectra of 
star-forming and E+A galaxies. Thus we can study the effect 
of age selections on stellar population synthesis. The results 
show that the stellar population synthesis does depend on the 
selections of ages. 

4. We fix the IMF, age, stellar evolution tracks, and change 
the selection of metallicities to construct more different SSP 
groups from the 6 EPS models. Then we fit the spectra of 
star-forming and E+A galaxies, in which way we can study 
the dependences of stellar population synthesis on the selec- 
tion of metallicities. The results show that the stellar popula- 
tion synthesis also depends on the selection of metallicities, 
but which is less important than ages. We notice that this 
less dependence on metallicity than age may be due to either 
the way we show our results or the classic age-metallicity 
degeneracy. 

5. We fix the age and metallicity and change the selection of 
stellar evolution tracks in BC03 and GALEV models respec- 
tively. Then we use them to fit the spectra of 6 types of galax- 
ies to check the effect of stellar evolution track on the stellar 
population synthesis. The results show that stellar evolution 
track does not affect much on the stellar population synthe- 
sis. 

6. We also compare the stellar population synthesis results 
among different types of galaxies, and suggest that there is a 
possible age sequence: the importance of young populations 
decreases from star-forming, composite, Seyfert 2, LINER 
to early-type galaxies, and E+A galaxies lie between com- 
posite galaxies and Seyfert 2s in most cases. 
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